3.1.17 \(\int (a+b x^2) (c+d x^2)^3 (e+f x^2)^2 \, dx\) [17]

3.1.17.1 Optimal result
3.1.17.2 Mathematica [A] (verified)
3.1.17.3 Rubi [A] (verified)
3.1.17.4 Maple [A] (verified)
3.1.17.5 Fricas [A] (verification not implemented)
3.1.17.6 Sympy [A] (verification not implemented)
3.1.17.7 Maxima [A] (verification not implemented)
3.1.17.8 Giac [A] (verification not implemented)
3.1.17.9 Mupad [B] (verification not implemented)

3.1.17.1 Optimal result

Integrand size = 26, antiderivative size = 226 \[ \int \left (a+b x^2\right ) \left (c+d x^2\right )^3 \left (e+f x^2\right )^2 \, dx=a c^3 e^2 x+\frac {1}{3} c^2 e (b c e+3 a d e+2 a c f) x^3+\frac {1}{5} c \left (b c e (3 d e+2 c f)+a \left (3 d^2 e^2+6 c d e f+c^2 f^2\right )\right ) x^5+\frac {1}{7} \left (b c \left (3 d^2 e^2+6 c d e f+c^2 f^2\right )+a d \left (d^2 e^2+6 c d e f+3 c^2 f^2\right )\right ) x^7+\frac {1}{9} d \left (a d f (2 d e+3 c f)+b \left (d^2 e^2+6 c d e f+3 c^2 f^2\right )\right ) x^9+\frac {1}{11} d^2 f (2 b d e+3 b c f+a d f) x^{11}+\frac {1}{13} b d^3 f^2 x^{13} \]

output
a*c^3*e^2*x+1/3*c^2*e*(2*a*c*f+3*a*d*e+b*c*e)*x^3+1/5*c*(b*c*e*(2*c*f+3*d* 
e)+a*(c^2*f^2+6*c*d*e*f+3*d^2*e^2))*x^5+1/7*(b*c*(c^2*f^2+6*c*d*e*f+3*d^2* 
e^2)+a*d*(3*c^2*f^2+6*c*d*e*f+d^2*e^2))*x^7+1/9*d*(a*d*f*(3*c*f+2*d*e)+b*( 
3*c^2*f^2+6*c*d*e*f+d^2*e^2))*x^9+1/11*d^2*f*(a*d*f+3*b*c*f+2*b*d*e)*x^11+ 
1/13*b*d^3*f^2*x^13
 
3.1.17.2 Mathematica [A] (verified)

Time = 0.06 (sec) , antiderivative size = 226, normalized size of antiderivative = 1.00 \[ \int \left (a+b x^2\right ) \left (c+d x^2\right )^3 \left (e+f x^2\right )^2 \, dx=a c^3 e^2 x+\frac {1}{3} c^2 e (b c e+3 a d e+2 a c f) x^3+\frac {1}{5} c \left (b c e (3 d e+2 c f)+a \left (3 d^2 e^2+6 c d e f+c^2 f^2\right )\right ) x^5+\frac {1}{7} \left (b c \left (3 d^2 e^2+6 c d e f+c^2 f^2\right )+a d \left (d^2 e^2+6 c d e f+3 c^2 f^2\right )\right ) x^7+\frac {1}{9} d \left (a d f (2 d e+3 c f)+b \left (d^2 e^2+6 c d e f+3 c^2 f^2\right )\right ) x^9+\frac {1}{11} d^2 f (2 b d e+3 b c f+a d f) x^{11}+\frac {1}{13} b d^3 f^2 x^{13} \]

input
Integrate[(a + b*x^2)*(c + d*x^2)^3*(e + f*x^2)^2,x]
 
output
a*c^3*e^2*x + (c^2*e*(b*c*e + 3*a*d*e + 2*a*c*f)*x^3)/3 + (c*(b*c*e*(3*d*e 
 + 2*c*f) + a*(3*d^2*e^2 + 6*c*d*e*f + c^2*f^2))*x^5)/5 + ((b*c*(3*d^2*e^2 
 + 6*c*d*e*f + c^2*f^2) + a*d*(d^2*e^2 + 6*c*d*e*f + 3*c^2*f^2))*x^7)/7 + 
(d*(a*d*f*(2*d*e + 3*c*f) + b*(d^2*e^2 + 6*c*d*e*f + 3*c^2*f^2))*x^9)/9 + 
(d^2*f*(2*b*d*e + 3*b*c*f + a*d*f)*x^11)/11 + (b*d^3*f^2*x^13)/13
 
3.1.17.3 Rubi [A] (verified)

Time = 0.44 (sec) , antiderivative size = 226, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.077, Rules used = {396, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \left (a+b x^2\right ) \left (c+d x^2\right )^3 \left (e+f x^2\right )^2 \, dx\)

\(\Big \downarrow \) 396

\(\displaystyle \int \left (d x^8 \left (a d f (3 c f+2 d e)+b \left (3 c^2 f^2+6 c d e f+d^2 e^2\right )\right )+x^6 \left (a d \left (3 c^2 f^2+6 c d e f+d^2 e^2\right )+b c \left (c^2 f^2+6 c d e f+3 d^2 e^2\right )\right )+c x^4 \left (a \left (c^2 f^2+6 c d e f+3 d^2 e^2\right )+b c e (2 c f+3 d e)\right )+c^2 e x^2 (2 a c f+3 a d e+b c e)+d^2 f x^{10} (a d f+3 b c f+2 b d e)+a c^3 e^2+b d^3 f^2 x^{12}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {1}{9} d x^9 \left (a d f (3 c f+2 d e)+b \left (3 c^2 f^2+6 c d e f+d^2 e^2\right )\right )+\frac {1}{7} x^7 \left (a d \left (3 c^2 f^2+6 c d e f+d^2 e^2\right )+b c \left (c^2 f^2+6 c d e f+3 d^2 e^2\right )\right )+\frac {1}{5} c x^5 \left (a \left (c^2 f^2+6 c d e f+3 d^2 e^2\right )+b c e (2 c f+3 d e)\right )+\frac {1}{3} c^2 e x^3 (2 a c f+3 a d e+b c e)+\frac {1}{11} d^2 f x^{11} (a d f+3 b c f+2 b d e)+a c^3 e^2 x+\frac {1}{13} b d^3 f^2 x^{13}\)

input
Int[(a + b*x^2)*(c + d*x^2)^3*(e + f*x^2)^2,x]
 
output
a*c^3*e^2*x + (c^2*e*(b*c*e + 3*a*d*e + 2*a*c*f)*x^3)/3 + (c*(b*c*e*(3*d*e 
 + 2*c*f) + a*(3*d^2*e^2 + 6*c*d*e*f + c^2*f^2))*x^5)/5 + ((b*c*(3*d^2*e^2 
 + 6*c*d*e*f + c^2*f^2) + a*d*(d^2*e^2 + 6*c*d*e*f + 3*c^2*f^2))*x^7)/7 + 
(d*(a*d*f*(2*d*e + 3*c*f) + b*(d^2*e^2 + 6*c*d*e*f + 3*c^2*f^2))*x^9)/9 + 
(d^2*f*(2*b*d*e + 3*b*c*f + a*d*f)*x^11)/11 + (b*d^3*f^2*x^13)/13
 

3.1.17.3.1 Defintions of rubi rules used

rule 396
Int[((a_) + (b_.)*(x_)^2)^(p_.)*((c_) + (d_.)*(x_)^2)^(q_.)*((e_) + (f_.)*( 
x_)^2)^(r_.), x_Symbol] :> Int[ExpandIntegrand[(a + b*x^2)^p*(c + d*x^2)^q* 
(e + f*x^2)^r, x], x] /; FreeQ[{a, b, c, d, e, f}, x] && IGtQ[p, 0] && IGtQ 
[q, 0] && IGtQ[r, 0]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.1.17.4 Maple [A] (verified)

Time = 3.28 (sec) , antiderivative size = 244, normalized size of antiderivative = 1.08

method result size
default \(\frac {b \,d^{3} f^{2} x^{13}}{13}+\frac {\left (\left (a \,d^{3}+3 b c \,d^{2}\right ) f^{2}+2 b \,d^{3} e f \right ) x^{11}}{11}+\frac {\left (\left (3 a c \,d^{2}+3 b \,c^{2} d \right ) f^{2}+2 \left (a \,d^{3}+3 b c \,d^{2}\right ) e f +b \,d^{3} e^{2}\right ) x^{9}}{9}+\frac {\left (\left (3 a \,c^{2} d +c^{3} b \right ) f^{2}+2 \left (3 a c \,d^{2}+3 b \,c^{2} d \right ) e f +\left (a \,d^{3}+3 b c \,d^{2}\right ) e^{2}\right ) x^{7}}{7}+\frac {\left (c^{3} a \,f^{2}+2 \left (3 a \,c^{2} d +c^{3} b \right ) e f +\left (3 a c \,d^{2}+3 b \,c^{2} d \right ) e^{2}\right ) x^{5}}{5}+\frac {\left (2 c^{3} a e f +\left (3 a \,c^{2} d +c^{3} b \right ) e^{2}\right ) x^{3}}{3}+a \,c^{3} e^{2} x\) \(244\)
norman \(\frac {b \,d^{3} f^{2} x^{13}}{13}+\left (\frac {1}{11} a \,d^{3} f^{2}+\frac {3}{11} b c \,d^{2} f^{2}+\frac {2}{11} b \,d^{3} e f \right ) x^{11}+\left (\frac {1}{3} a c \,d^{2} f^{2}+\frac {2}{9} a \,d^{3} e f +\frac {1}{3} b \,c^{2} d \,f^{2}+\frac {2}{3} b c \,d^{2} e f +\frac {1}{9} b \,d^{3} e^{2}\right ) x^{9}+\left (\frac {3}{7} a \,c^{2} d \,f^{2}+\frac {6}{7} a c \,d^{2} e f +\frac {1}{7} a \,d^{3} e^{2}+\frac {1}{7} b \,c^{3} f^{2}+\frac {6}{7} b \,c^{2} d e f +\frac {3}{7} b c \,d^{2} e^{2}\right ) x^{7}+\left (\frac {1}{5} c^{3} a \,f^{2}+\frac {6}{5} a \,c^{2} d e f +\frac {3}{5} a c \,d^{2} e^{2}+\frac {2}{5} b \,c^{3} e f +\frac {3}{5} b \,c^{2} d \,e^{2}\right ) x^{5}+\left (\frac {2}{3} c^{3} a e f +a \,c^{2} d \,e^{2}+\frac {1}{3} b \,c^{3} e^{2}\right ) x^{3}+a \,c^{3} e^{2} x\) \(249\)
gosper \(\frac {1}{13} b \,d^{3} f^{2} x^{13}+\frac {1}{11} x^{11} a \,d^{3} f^{2}+\frac {3}{11} x^{11} b c \,d^{2} f^{2}+\frac {2}{11} x^{11} b \,d^{3} e f +\frac {1}{3} x^{9} a c \,d^{2} f^{2}+\frac {2}{9} x^{9} a \,d^{3} e f +\frac {1}{3} x^{9} b \,c^{2} d \,f^{2}+\frac {2}{3} x^{9} b c \,d^{2} e f +\frac {1}{9} x^{9} b \,d^{3} e^{2}+\frac {3}{7} x^{7} a \,c^{2} d \,f^{2}+\frac {6}{7} x^{7} a c \,d^{2} e f +\frac {1}{7} x^{7} a \,d^{3} e^{2}+\frac {1}{7} x^{7} b \,c^{3} f^{2}+\frac {6}{7} x^{7} b \,c^{2} d e f +\frac {3}{7} x^{7} b c \,d^{2} e^{2}+\frac {1}{5} x^{5} c^{3} a \,f^{2}+\frac {6}{5} x^{5} a \,c^{2} d e f +\frac {3}{5} x^{5} a c \,d^{2} e^{2}+\frac {2}{5} x^{5} b \,c^{3} e f +\frac {3}{5} x^{5} b \,c^{2} d \,e^{2}+\frac {2}{3} x^{3} c^{3} a e f +x^{3} a \,c^{2} d \,e^{2}+\frac {1}{3} x^{3} b \,c^{3} e^{2}+a \,c^{3} e^{2} x\) \(290\)
risch \(\frac {1}{13} b \,d^{3} f^{2} x^{13}+\frac {1}{11} x^{11} a \,d^{3} f^{2}+\frac {3}{11} x^{11} b c \,d^{2} f^{2}+\frac {2}{11} x^{11} b \,d^{3} e f +\frac {1}{3} x^{9} a c \,d^{2} f^{2}+\frac {2}{9} x^{9} a \,d^{3} e f +\frac {1}{3} x^{9} b \,c^{2} d \,f^{2}+\frac {2}{3} x^{9} b c \,d^{2} e f +\frac {1}{9} x^{9} b \,d^{3} e^{2}+\frac {3}{7} x^{7} a \,c^{2} d \,f^{2}+\frac {6}{7} x^{7} a c \,d^{2} e f +\frac {1}{7} x^{7} a \,d^{3} e^{2}+\frac {1}{7} x^{7} b \,c^{3} f^{2}+\frac {6}{7} x^{7} b \,c^{2} d e f +\frac {3}{7} x^{7} b c \,d^{2} e^{2}+\frac {1}{5} x^{5} c^{3} a \,f^{2}+\frac {6}{5} x^{5} a \,c^{2} d e f +\frac {3}{5} x^{5} a c \,d^{2} e^{2}+\frac {2}{5} x^{5} b \,c^{3} e f +\frac {3}{5} x^{5} b \,c^{2} d \,e^{2}+\frac {2}{3} x^{3} c^{3} a e f +x^{3} a \,c^{2} d \,e^{2}+\frac {1}{3} x^{3} b \,c^{3} e^{2}+a \,c^{3} e^{2} x\) \(290\)
parallelrisch \(\frac {1}{13} b \,d^{3} f^{2} x^{13}+\frac {1}{11} x^{11} a \,d^{3} f^{2}+\frac {3}{11} x^{11} b c \,d^{2} f^{2}+\frac {2}{11} x^{11} b \,d^{3} e f +\frac {1}{3} x^{9} a c \,d^{2} f^{2}+\frac {2}{9} x^{9} a \,d^{3} e f +\frac {1}{3} x^{9} b \,c^{2} d \,f^{2}+\frac {2}{3} x^{9} b c \,d^{2} e f +\frac {1}{9} x^{9} b \,d^{3} e^{2}+\frac {3}{7} x^{7} a \,c^{2} d \,f^{2}+\frac {6}{7} x^{7} a c \,d^{2} e f +\frac {1}{7} x^{7} a \,d^{3} e^{2}+\frac {1}{7} x^{7} b \,c^{3} f^{2}+\frac {6}{7} x^{7} b \,c^{2} d e f +\frac {3}{7} x^{7} b c \,d^{2} e^{2}+\frac {1}{5} x^{5} c^{3} a \,f^{2}+\frac {6}{5} x^{5} a \,c^{2} d e f +\frac {3}{5} x^{5} a c \,d^{2} e^{2}+\frac {2}{5} x^{5} b \,c^{3} e f +\frac {3}{5} x^{5} b \,c^{2} d \,e^{2}+\frac {2}{3} x^{3} c^{3} a e f +x^{3} a \,c^{2} d \,e^{2}+\frac {1}{3} x^{3} b \,c^{3} e^{2}+a \,c^{3} e^{2} x\) \(290\)

input
int((b*x^2+a)*(d*x^2+c)^3*(f*x^2+e)^2,x,method=_RETURNVERBOSE)
 
output
1/13*b*d^3*f^2*x^13+1/11*((a*d^3+3*b*c*d^2)*f^2+2*b*d^3*e*f)*x^11+1/9*((3* 
a*c*d^2+3*b*c^2*d)*f^2+2*(a*d^3+3*b*c*d^2)*e*f+b*d^3*e^2)*x^9+1/7*((3*a*c^ 
2*d+b*c^3)*f^2+2*(3*a*c*d^2+3*b*c^2*d)*e*f+(a*d^3+3*b*c*d^2)*e^2)*x^7+1/5* 
(c^3*a*f^2+2*(3*a*c^2*d+b*c^3)*e*f+(3*a*c*d^2+3*b*c^2*d)*e^2)*x^5+1/3*(2*c 
^3*a*e*f+(3*a*c^2*d+b*c^3)*e^2)*x^3+a*c^3*e^2*x
 
3.1.17.5 Fricas [A] (verification not implemented)

Time = 0.33 (sec) , antiderivative size = 239, normalized size of antiderivative = 1.06 \[ \int \left (a+b x^2\right ) \left (c+d x^2\right )^3 \left (e+f x^2\right )^2 \, dx=\frac {1}{13} \, b d^{3} f^{2} x^{13} + \frac {1}{11} \, {\left (2 \, b d^{3} e f + {\left (3 \, b c d^{2} + a d^{3}\right )} f^{2}\right )} x^{11} + \frac {1}{9} \, {\left (b d^{3} e^{2} + 2 \, {\left (3 \, b c d^{2} + a d^{3}\right )} e f + 3 \, {\left (b c^{2} d + a c d^{2}\right )} f^{2}\right )} x^{9} + \frac {1}{7} \, {\left ({\left (3 \, b c d^{2} + a d^{3}\right )} e^{2} + 6 \, {\left (b c^{2} d + a c d^{2}\right )} e f + {\left (b c^{3} + 3 \, a c^{2} d\right )} f^{2}\right )} x^{7} + a c^{3} e^{2} x + \frac {1}{5} \, {\left (a c^{3} f^{2} + 3 \, {\left (b c^{2} d + a c d^{2}\right )} e^{2} + 2 \, {\left (b c^{3} + 3 \, a c^{2} d\right )} e f\right )} x^{5} + \frac {1}{3} \, {\left (2 \, a c^{3} e f + {\left (b c^{3} + 3 \, a c^{2} d\right )} e^{2}\right )} x^{3} \]

input
integrate((b*x^2+a)*(d*x^2+c)^3*(f*x^2+e)^2,x, algorithm="fricas")
 
output
1/13*b*d^3*f^2*x^13 + 1/11*(2*b*d^3*e*f + (3*b*c*d^2 + a*d^3)*f^2)*x^11 + 
1/9*(b*d^3*e^2 + 2*(3*b*c*d^2 + a*d^3)*e*f + 3*(b*c^2*d + a*c*d^2)*f^2)*x^ 
9 + 1/7*((3*b*c*d^2 + a*d^3)*e^2 + 6*(b*c^2*d + a*c*d^2)*e*f + (b*c^3 + 3* 
a*c^2*d)*f^2)*x^7 + a*c^3*e^2*x + 1/5*(a*c^3*f^2 + 3*(b*c^2*d + a*c*d^2)*e 
^2 + 2*(b*c^3 + 3*a*c^2*d)*e*f)*x^5 + 1/3*(2*a*c^3*e*f + (b*c^3 + 3*a*c^2* 
d)*e^2)*x^3
 
3.1.17.6 Sympy [A] (verification not implemented)

Time = 0.04 (sec) , antiderivative size = 304, normalized size of antiderivative = 1.35 \[ \int \left (a+b x^2\right ) \left (c+d x^2\right )^3 \left (e+f x^2\right )^2 \, dx=a c^{3} e^{2} x + \frac {b d^{3} f^{2} x^{13}}{13} + x^{11} \left (\frac {a d^{3} f^{2}}{11} + \frac {3 b c d^{2} f^{2}}{11} + \frac {2 b d^{3} e f}{11}\right ) + x^{9} \left (\frac {a c d^{2} f^{2}}{3} + \frac {2 a d^{3} e f}{9} + \frac {b c^{2} d f^{2}}{3} + \frac {2 b c d^{2} e f}{3} + \frac {b d^{3} e^{2}}{9}\right ) + x^{7} \cdot \left (\frac {3 a c^{2} d f^{2}}{7} + \frac {6 a c d^{2} e f}{7} + \frac {a d^{3} e^{2}}{7} + \frac {b c^{3} f^{2}}{7} + \frac {6 b c^{2} d e f}{7} + \frac {3 b c d^{2} e^{2}}{7}\right ) + x^{5} \left (\frac {a c^{3} f^{2}}{5} + \frac {6 a c^{2} d e f}{5} + \frac {3 a c d^{2} e^{2}}{5} + \frac {2 b c^{3} e f}{5} + \frac {3 b c^{2} d e^{2}}{5}\right ) + x^{3} \cdot \left (\frac {2 a c^{3} e f}{3} + a c^{2} d e^{2} + \frac {b c^{3} e^{2}}{3}\right ) \]

input
integrate((b*x**2+a)*(d*x**2+c)**3*(f*x**2+e)**2,x)
 
output
a*c**3*e**2*x + b*d**3*f**2*x**13/13 + x**11*(a*d**3*f**2/11 + 3*b*c*d**2* 
f**2/11 + 2*b*d**3*e*f/11) + x**9*(a*c*d**2*f**2/3 + 2*a*d**3*e*f/9 + b*c* 
*2*d*f**2/3 + 2*b*c*d**2*e*f/3 + b*d**3*e**2/9) + x**7*(3*a*c**2*d*f**2/7 
+ 6*a*c*d**2*e*f/7 + a*d**3*e**2/7 + b*c**3*f**2/7 + 6*b*c**2*d*e*f/7 + 3* 
b*c*d**2*e**2/7) + x**5*(a*c**3*f**2/5 + 6*a*c**2*d*e*f/5 + 3*a*c*d**2*e** 
2/5 + 2*b*c**3*e*f/5 + 3*b*c**2*d*e**2/5) + x**3*(2*a*c**3*e*f/3 + a*c**2* 
d*e**2 + b*c**3*e**2/3)
 
3.1.17.7 Maxima [A] (verification not implemented)

Time = 0.22 (sec) , antiderivative size = 239, normalized size of antiderivative = 1.06 \[ \int \left (a+b x^2\right ) \left (c+d x^2\right )^3 \left (e+f x^2\right )^2 \, dx=\frac {1}{13} \, b d^{3} f^{2} x^{13} + \frac {1}{11} \, {\left (2 \, b d^{3} e f + {\left (3 \, b c d^{2} + a d^{3}\right )} f^{2}\right )} x^{11} + \frac {1}{9} \, {\left (b d^{3} e^{2} + 2 \, {\left (3 \, b c d^{2} + a d^{3}\right )} e f + 3 \, {\left (b c^{2} d + a c d^{2}\right )} f^{2}\right )} x^{9} + \frac {1}{7} \, {\left ({\left (3 \, b c d^{2} + a d^{3}\right )} e^{2} + 6 \, {\left (b c^{2} d + a c d^{2}\right )} e f + {\left (b c^{3} + 3 \, a c^{2} d\right )} f^{2}\right )} x^{7} + a c^{3} e^{2} x + \frac {1}{5} \, {\left (a c^{3} f^{2} + 3 \, {\left (b c^{2} d + a c d^{2}\right )} e^{2} + 2 \, {\left (b c^{3} + 3 \, a c^{2} d\right )} e f\right )} x^{5} + \frac {1}{3} \, {\left (2 \, a c^{3} e f + {\left (b c^{3} + 3 \, a c^{2} d\right )} e^{2}\right )} x^{3} \]

input
integrate((b*x^2+a)*(d*x^2+c)^3*(f*x^2+e)^2,x, algorithm="maxima")
 
output
1/13*b*d^3*f^2*x^13 + 1/11*(2*b*d^3*e*f + (3*b*c*d^2 + a*d^3)*f^2)*x^11 + 
1/9*(b*d^3*e^2 + 2*(3*b*c*d^2 + a*d^3)*e*f + 3*(b*c^2*d + a*c*d^2)*f^2)*x^ 
9 + 1/7*((3*b*c*d^2 + a*d^3)*e^2 + 6*(b*c^2*d + a*c*d^2)*e*f + (b*c^3 + 3* 
a*c^2*d)*f^2)*x^7 + a*c^3*e^2*x + 1/5*(a*c^3*f^2 + 3*(b*c^2*d + a*c*d^2)*e 
^2 + 2*(b*c^3 + 3*a*c^2*d)*e*f)*x^5 + 1/3*(2*a*c^3*e*f + (b*c^3 + 3*a*c^2* 
d)*e^2)*x^3
 
3.1.17.8 Giac [A] (verification not implemented)

Time = 0.30 (sec) , antiderivative size = 289, normalized size of antiderivative = 1.28 \[ \int \left (a+b x^2\right ) \left (c+d x^2\right )^3 \left (e+f x^2\right )^2 \, dx=\frac {1}{13} \, b d^{3} f^{2} x^{13} + \frac {2}{11} \, b d^{3} e f x^{11} + \frac {3}{11} \, b c d^{2} f^{2} x^{11} + \frac {1}{11} \, a d^{3} f^{2} x^{11} + \frac {1}{9} \, b d^{3} e^{2} x^{9} + \frac {2}{3} \, b c d^{2} e f x^{9} + \frac {2}{9} \, a d^{3} e f x^{9} + \frac {1}{3} \, b c^{2} d f^{2} x^{9} + \frac {1}{3} \, a c d^{2} f^{2} x^{9} + \frac {3}{7} \, b c d^{2} e^{2} x^{7} + \frac {1}{7} \, a d^{3} e^{2} x^{7} + \frac {6}{7} \, b c^{2} d e f x^{7} + \frac {6}{7} \, a c d^{2} e f x^{7} + \frac {1}{7} \, b c^{3} f^{2} x^{7} + \frac {3}{7} \, a c^{2} d f^{2} x^{7} + \frac {3}{5} \, b c^{2} d e^{2} x^{5} + \frac {3}{5} \, a c d^{2} e^{2} x^{5} + \frac {2}{5} \, b c^{3} e f x^{5} + \frac {6}{5} \, a c^{2} d e f x^{5} + \frac {1}{5} \, a c^{3} f^{2} x^{5} + \frac {1}{3} \, b c^{3} e^{2} x^{3} + a c^{2} d e^{2} x^{3} + \frac {2}{3} \, a c^{3} e f x^{3} + a c^{3} e^{2} x \]

input
integrate((b*x^2+a)*(d*x^2+c)^3*(f*x^2+e)^2,x, algorithm="giac")
 
output
1/13*b*d^3*f^2*x^13 + 2/11*b*d^3*e*f*x^11 + 3/11*b*c*d^2*f^2*x^11 + 1/11*a 
*d^3*f^2*x^11 + 1/9*b*d^3*e^2*x^9 + 2/3*b*c*d^2*e*f*x^9 + 2/9*a*d^3*e*f*x^ 
9 + 1/3*b*c^2*d*f^2*x^9 + 1/3*a*c*d^2*f^2*x^9 + 3/7*b*c*d^2*e^2*x^7 + 1/7* 
a*d^3*e^2*x^7 + 6/7*b*c^2*d*e*f*x^7 + 6/7*a*c*d^2*e*f*x^7 + 1/7*b*c^3*f^2* 
x^7 + 3/7*a*c^2*d*f^2*x^7 + 3/5*b*c^2*d*e^2*x^5 + 3/5*a*c*d^2*e^2*x^5 + 2/ 
5*b*c^3*e*f*x^5 + 6/5*a*c^2*d*e*f*x^5 + 1/5*a*c^3*f^2*x^5 + 1/3*b*c^3*e^2* 
x^3 + a*c^2*d*e^2*x^3 + 2/3*a*c^3*e*f*x^3 + a*c^3*e^2*x
 
3.1.17.9 Mupad [B] (verification not implemented)

Time = 5.18 (sec) , antiderivative size = 233, normalized size of antiderivative = 1.03 \[ \int \left (a+b x^2\right ) \left (c+d x^2\right )^3 \left (e+f x^2\right )^2 \, dx=x^5\,\left (\frac {2\,b\,c^3\,e\,f}{5}+\frac {a\,c^3\,f^2}{5}+\frac {3\,b\,c^2\,d\,e^2}{5}+\frac {6\,a\,c^2\,d\,e\,f}{5}+\frac {3\,a\,c\,d^2\,e^2}{5}\right )+x^9\,\left (\frac {b\,c^2\,d\,f^2}{3}+\frac {2\,b\,c\,d^2\,e\,f}{3}+\frac {a\,c\,d^2\,f^2}{3}+\frac {b\,d^3\,e^2}{9}+\frac {2\,a\,d^3\,e\,f}{9}\right )+x^7\,\left (\frac {b\,c^3\,f^2}{7}+\frac {6\,b\,c^2\,d\,e\,f}{7}+\frac {3\,a\,c^2\,d\,f^2}{7}+\frac {3\,b\,c\,d^2\,e^2}{7}+\frac {6\,a\,c\,d^2\,e\,f}{7}+\frac {a\,d^3\,e^2}{7}\right )+\frac {b\,d^3\,f^2\,x^{13}}{13}+\frac {c^2\,e\,x^3\,\left (2\,a\,c\,f+3\,a\,d\,e+b\,c\,e\right )}{3}+\frac {d^2\,f\,x^{11}\,\left (a\,d\,f+3\,b\,c\,f+2\,b\,d\,e\right )}{11}+a\,c^3\,e^2\,x \]

input
int((a + b*x^2)*(c + d*x^2)^3*(e + f*x^2)^2,x)
 
output
x^5*((a*c^3*f^2)/5 + (2*b*c^3*e*f)/5 + (3*a*c*d^2*e^2)/5 + (3*b*c^2*d*e^2) 
/5 + (6*a*c^2*d*e*f)/5) + x^9*((b*d^3*e^2)/9 + (2*a*d^3*e*f)/9 + (a*c*d^2* 
f^2)/3 + (b*c^2*d*f^2)/3 + (2*b*c*d^2*e*f)/3) + x^7*((a*d^3*e^2)/7 + (b*c^ 
3*f^2)/7 + (3*a*c^2*d*f^2)/7 + (3*b*c*d^2*e^2)/7 + (6*a*c*d^2*e*f)/7 + (6* 
b*c^2*d*e*f)/7) + (b*d^3*f^2*x^13)/13 + (c^2*e*x^3*(2*a*c*f + 3*a*d*e + b* 
c*e))/3 + (d^2*f*x^11*(a*d*f + 3*b*c*f + 2*b*d*e))/11 + a*c^3*e^2*x